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Abstract: We analyze vacuum decay in brane world setups, where a free scalar field in five 
dimensions has a localized potential admitting metastable vacua. We study in particular 
the bounce solution and its properties in flat and warped spaces. In the latter case, placing 
into a deeply warped region the term in the potential that lifts the vacuum degeneracy, 
can increase indefinitely the lifetime of the false vacuum. We discuss the application to 
metastable vacua in supersymmetric brane-world constructions. 
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1. Introduction and conclusions 

The fate of metastable vacua in field theory [|l|, ^ is of great interest in cosmology and 
particle physics. The dynamics of their quantum decay toward the true vacuum rely on 
the knowledge of the classical Euclidean "bounce" solution. The study of finite energy 
soliton solutions in Minkowski space and of finite action solutions in Euclidean space is 
therefore crucial for the study of metastable vacua. The purpose of this paper is to study 
the generalization of this problem in the case of vacua generated by a scalar field living in 
a higher-dimensional spacetime, with a scalar potential localized in four dimensions. The 
field is therefore free in the bulk, with the scalar potential generating non-trivial boundary 
conditions. One of our main motivation for studying the case of a boundary potential is 
the generalization to supersymmetric theories with metastable vacua |^. In this case, 
constraints coming from higher-dimensional supersymmetry are such that it is much easier 
to construct models with localized (as opposed to bulk) superpotentials. 

The search for classical solutions in this case turns out to be very interesting and rich. 
A first question is the dependence of the classical solution on the geometry of the internal 
space, flat or warped, and on its size. In the case of the bounce, we would like to understand 
the dependence of the width of the wall in the thin wall approximation on the extra 
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dimension. Moreover, since all nontrivial dynamics is encoded in boundary conditions, it 
suggests that the problem could be tractable to some extent even in the case where the 
internal space is warped. In this last case, there are several interesting questions arising. 
First of all, even if we are in a regime in which the 4d effective theory is valid, i.e. there is a 
mode much lighter than the KK masses, it is possible that the barrier separating the false 
from the true minimum is much higher than the mass of the lowest-lying KK states. In 
this case, despite the validity of the 4d effective action describing the lightest mode, there 
is no a priori reason why the classical solution should be the standard 4d one. Secondly, 
it is reasonable to expect that by placing the term lifting the degeneracy between the true 
and the false vacuum into a deeply warped region, it will be redshifted to small values, thus 
increasing indefinitely the lifetime of the false vacuum. If this were indeed possible, there 
would be no practical difference between living in the true vacuum or the false vacuum ! 
In this paper, we will be able to answer some of these questions, whereas other questions 
will be addressed only partially and will need future work for a complete understanding. 

The plan of the paper is as follows. In Section ^ we discuss vacuum decay in a toy-model 
consisting of a single scalar field in a spacetime with one flat, compact extra-dimension. 
We show that when there is a light (compared to the KK scale) mode, the effective theory 
is precisely the one which admits a standard kink solution. We then work out the 5d analog 
of the classical 4d field equation. This can be written as a 4d differential equation, which 
allows a systematic calculation of the corrections to the 4d kink solution. The equation 
involves a differential operator containing higher derivative terms. We work out the size of 
the kink in the large radius limit, show that the kink become broader with increasing the 
radius and study its behavior near the origin. We use this to write the equation defining 
the euclidian bounce and check the validity of the thin wall approximation, finding that it 
gets worse in the 5d limit. 

In Section |^ we perform the same analysis including the effect of the warping of the 
extra-dimension. We find that the 4d limit and the thin-wall approximation become very 
accurate because of the warping. However Coleman-de Luccia gravitational effects can 
become important and even completely lock the decay of the false vacuum. 

In Section ^ we extend these considerations to the supersymmetric case. Motivated by 
the D3/D7 brane realizations of the ISS model, we discuss the AdS^ version of the ISS 
model HI, with ISS gauge group and quarks living on the UV boundary and the mesons 
living in the 5d bulk. The meson-quark coupling is then localized on the UV brane, whereas 
the mesonic linear term in superpotential is put on the IR boundary. We find that, due 
to the warping, the (mass)^ parameter is naturally redshifted to small values, whereas 
metastable supersymmetry breaking becomes a non-local (in the extra dimension) effect. 
This has again the net effect of increasing correspondingly the lifetime of the metastable 
vacuum. We also analyze briefly the case where the whole superpotential is localized on 
the UV brane and a light mode is achieved by adding a bulk mass term for the mesons 
hypermultiplet. In this case mass scales are redshifted again due to a different effect, the 
value of the mesonic wave function on the UV brane. Both examples have a natural 4d 
holographic interpretation via the AdS/CFT correspondence. 

Some technical details of the computations are left to three Appendices. 
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2. Flat 4+1 Dimensions 



In this section we consider a massless scalar field in a 4+1 dimensional flat spacetime^ in 
which the 5th direction (labeled by the coordinate y) extends between two rigid branes at 
y = and y = ttR. The bulk action is that of a free massless field, all nontrivial potential 
terms appearing on the boundaries'^ : 

S = -^l d^xdydA^d^^- j d^xVom\y=, + J d^xV^my^^^. (2.1) 
The field equations and boundary conditions read: 

d^^ + a^a^«> = 0, (2.2) 

^y%^o = ^^ (2-3) 

5 $1 = (2.4) 

2.1 The Kink 

Consider the situation where the brane potentials are given by: 

Vom = \{^^-v^)\ V,{^) = Q (2.5) 

From eqs. ( |2.2[j2.4| ) we see immediately that there are two "vacuum" solutions ^±{x^y) = 
^v. One can ask whether there exist a solution interpolating between the two vacua, 
analogous to the four-dimensional domain wall (kink) that one finds with the same quartic 
potential (see Appendix^). 

Notice that, since $ is canonically normalized in 5D, and has mass dimension 3/2, the 
parameters in ( |2.5D have unusual mass dimensions: 

[A] = M-2, [v] = M^/\ (2.6) 

2.1.1 Effective 4D theory 

A kink-like solution is expected to exist at least in a certain region of parameter space, 
where one can give a four-dimensional effective description of the model. To see this, 
consider the linearized fluctuations around one of the two vacua (say $_): 

^{x,y) = -V + 6^{x,y). (2.7) 

Decomposing the solution in eigenstates of the 4D D'Alambertian, 5^{x,y) = 4>{y)xix), 
^4:X{y) = f^'^x{v)i the mass spectrum is obtained by linearizing the boundary conditions 

[^y^ = ^^l^]y=,, [dy^ = %=^n (2.8) 



^We use signature ( — h + + +)• Throughout the paper we consistently neglect the backreaction of the 
scalar field on the geometry ^ 

^The case of the bulk potential and the kink solution in the extra coordinate did lead historically to the 
first brane world proposals [M. For bounce solutions for brane localized fields, see e.g. fel. 
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where /ig = 2Au^. The mass eigenstates are the solutions of the equation: 

mtanrmrR = /Xq (2-9) 

and the profile wave-function for a given mode of mass m is: 

(j)m{y) = cos[m{y - ttR)]. (2.10) 

We have a low-energy, 4D effective theory for the lowest-lying mode (of mass mg) if rnoR ^ 
1. This description is valid for energies much smaller than the mass of the next KK mode, 
which is of order 1/R. Under this conditions we can expand the tangent in eq. (|2.9| ) and 
obtain: 

and the condition for the existence of a 4D description reads, in terms of the original 
parameters of the model: 

2RXv'^ < vr . (2.12) 

Under these condition, inserting ^(x,y) = —v + 4>oiy)xoix) in the original action and 
integrating over y, we obtain the low-energy 4D effective action for the lowest-lying mode 
Xo{x)- After some integration by parts and using the bulk field equation we obtain: 



Seff = I d X 
where the effective potential is: 



(2.13) 



2 9 3,^. 



Veff = -fxt - |X8 + jxh (2.14) 



,2 



2Xv^ 3Xv 



h = ^-^. (2.15) 



°^ vri? ' ^ (7ri2)3/2' (7ri?)2- 

The extra factors of vri? in the effective parameters come from the normalized wave- function 
profile 0o(y) = (l/VvrR) cos mo(y — vri?) evaluated in ?/ = 0. It is easy to check that 
the potential ( 2.14D has two zero-energy minima at xo = 0, 2wV ttR and a maximum at 
Xo = vy/ ttR with V{vV ttR) = /A. In terms of the original field $ = —v + 4>oXf) these 
correspond exactly to the original two minima at $ = ztv and maximum at $ = 0. 

Due to the standard double-well form of the effective potential the field equation derived 
from the effective action ( 2.13| ), 

dlxix) = ^ (2.16) 



admits a kink solutions interpolating between the two vacua Xo = and xo = 2^%/ R, which 
according to eqs. ( [A.2| , [A. 3D has the form: 

Xkink{x) = uV7rR(l -Ftanh/ia;) , A*^ = 2 (^^VTrRj = -^^^ = (2.17) 

The kink energy density is of the order of the height of the potential barrier, Af^. In 
order for the solution we found to be reliable, this energy density must be below the KK 
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scale, we thus have the additional requirement Xv^ <^ l/R^. This, together with ( 2.12| ), 
sets the range of validity of the kink solution we found. A sufficient condition is : 



V < R-^/^, X < R^. (2.18) 

Although the energy density of the domain wall is i2-independent, the integrated total 
energy is not: 

E = 1^"" dx(^^^y r^v'7rRf, = Vx^v^ (2.19) 
2.1.2 The 5D equation 

The argument of the previous subsection suggest that a kink solution to the model (I 



should exists, at least in the range of parameters satisfying ( 2.18 ). Now we want to look for 
similar solutions from a purely 5D perspective, without having to rely on the 4D effective 
theory approach. 

Let us return to eqs. ( p.2| - |2.4D . We look for solutions depending on y and one of the 
Minkowski coordinates (say x). The most general (real) solution to ( |2.2D can be written 
as: 

<^{x, y) = g{x + iy) + {g{x + iy))* . (2.20) 
The boundary condition at y = ttR tells us that 

Im [g'{x + mR)\ = 0, (2.21) 

where a prime denotes derivative w.r.t. the argument. This equation is satisfied if F[x) = 
g{x + iirR) is a real function'^. This also implies that, for z complex, {F{z))* = F(z*). 

Next, consider the boundary conditions at ?/ = 0. Defining f{x) = ^{x,0), h{x) = 
9y$(x,0), eq. dU) reads: 

h{x) = Xfix) [if{x)f - v^] . (2.22) 

Formally, we can write: 

fix) = g{x) + g{x)* = F{x - iirR) + F{x + iirR) 

= (^expl-iTTRda:] + exp[iTrRd^]jF{x), (2.23) 
h{x) = ig{x) — ig {x)* = iF'{x — iirR) — iF'{x + iirR) 

= i(^exp[—i'irRdx] — expliirRdxijdxFix) 

= tan{TTRdx)dxf{x). (2.24) 



^That is, g{z) has an expansion of the form 
with real coefficients c„. 
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Using the last hne in eq. ( ^.24 ) we arrive at a closed equation^ for f{x): 



(2.26) 



Another, maybe less general way, in order to arrive at ( 2.26| ) is to start from the bulk 
solution 

^{x,y) = j dpap eP"" cos{py + ap) , (2.27) 

where Up (ap) are arbitrary coefficients (phases). Boundary conditions at y = ttR fixes 
ap = —p ttR, whereas boundary conditions in y = gives by a straightforward computation 



(2.26), by using the replacement p ^ dx- This method will generalize in a straightforward 
manner to the warped case discussed in the next section. 

A solution to ( 2.26D gives f{x) = <^(x, 0), which can then be extended into the bulk to a 
full solution: 



<I>(x, y) = Re[f{x + iy)] + tax).{TTRdx)Im[f {x + iy)] 
= [cosydx + {iaxi-KRdx){sinydx)\ f{x) . 



(2.28) 



Eq. (|2]2|) has various interesting properties. It should be understood as a series in 



derivatives of increasing order. If we take f{x) to be a 4D mass eigenstate, f{x) 



and linearize the r.h.s, we get back to the eigenvalue equation ( |2.9| ). So the information 
about the spectrum of the model is contained in ( |2.26| ). 

Now suppose that we can keep the lowest order in the expansion of the l.h.s. (for any 
given solution we can later check whether this approximation is justified). We get a second 
order equation for / which looks exactly as the one for the 4D kink: 



A 



(JxJ ^J^J [J 



whose solution is again given by eq. ([A. 3]) : 

f{x) = V tanh^x. 



/i^ = At;V(2vri?). 



(2.29) 



(2.30) 



Notice that the characteristic scale fi is the same as in eq. (2.17). 



We can extract considerable information from eq. ( 2.26 ) even when the 4D limit does 
not hold. Consider a solution /(x) that approaches if as x — > ±00. We can estimate the 
width of the kink by expanding /(x) = —v + ry(x) and solving the linear equation for rj 



^As for the standard kink, this equation can be obtained from a the point-particle analog model, with 
potential —V{f) and an exotic kinetic term; 

1 



dx 



V{.f)- - ftan{7vRd^)d^f 



(2.25) 



which is the same as the "effective action" whose variation gives (2.26) 
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in the asymptotic large \x\ region: assuming rj{x) ~ e l^l/'*", where Iw is a measure of the 
wall width, we find: 

1 /vri?\ „. 2 



tan — = 2\v^ . (2.3f ) 

'"W \ '"W J 

In the 4D limit we get the expected result, namely 1^ = 1//^, with /i as in ( p. 30 ). In any 
case, the size of the wall cannot exceed 2i?, and this value is approached in the opposite 
limit, when R\v^ >> 1. 

Another interesting length scale is the one corresponding to the regime of the validity 
of the linear slope of the solution f{x) ~ [1/Iq)x^ when the variation (derivative) of the 
field / is maximal^. This can be estimated by linearizing eq. ( |2.26D around / = 0. Setting 
f{x) = 7/sin(x//o), with rj a constant, we find: 

- tanh — = Xv^ . (2.32) 



^0 V ^0 

In the 4d limit we get as expected /q ~ I/a*) whereas in the 5d regime we get /q ~ (l/Au^). 
In the 5d limit therefore, the size of the kink becomes larger and larger, whereas large 
variations of the field are confined into a fixed region. 
We can put eq. ( |2.26| ) in integral form. Going in Fourier space, and using the identity 

and /(O) = 0, we obtain: 

/(n) = ^jdt log [sinhl^^l] f{t){f{tf - 1) (2.34) 

where we have defined f{u) = f{u/{\v'^)). Since the kernel in ( p.34| ) behaves as \x — 1\ 
for large t, f{t) must necessarily approach one of the extrema / = 0, ztw as t — > ±00. 

It is also of interest to find the generalization of equ. ( [A.4p which expresses the vanishing 
of the "energy" of the point particle analog system. The determination of this quantity 
turns out to be non-trivial. The details of the calculation are relegated to the Appendix 
where we show that the kink verifies the following equation 

£ = ^-^^^ m? - {tM^Rd)dff] - -^Vif) = 0. (2.35) 
The leading terms in the expansion in powers of R are 

kdff-^V{f) + ••■ = (). (2.36) 



They reproduce the 4D equation ( [A.4| ). The first order correction to the 4D equation of 
motion can also be simply obtained: 

+ l{7rRfi2d'fdf - {d'ff)] - -^VU) = 0. (2.37) 



In the usual d + 1 kink solution (A.3) both 1^ and lo are of order 1/fi. 
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If we expand / in powers of R and write f = fo + {t^R)"^ f2 + • • • , then /2 can be determined 
from the first order equation: 



/0/2 



0, 



(2.38) 



where /o is the 4D solution viax\h.^x. Using the zeroth order equation of motion /q 
;^y(/o) the solution can be written as 



f2{x) 



ax) 



6 



{dkf 



(2.39) 



where the integration constant was fixed by requiring /2 to be odd. Finally the integration 
can be done to yield 



1 



h{x) = -ax) 



6" 



-^.^x + 4^-/o(x) 



(2.40) 



Explicitly, the solution to the first nontrivial order reads 

^2 



f{x) = f tanh(/ux) 



^ (27ri?/i)^ / /ix 
6cosh^(/ia;) Vtanh(^x) 



(2.41) 



This shows that the zeroth order approximation is valid as long as (27ri?/i)^ is much smaller 
than one. 

The identification of the kink as a solution of the equation = 0, has an important 
consequence: such solutions can never cross the lines / = ztv. This is analog to the usual 
two-derivative kink: there, eq. ( [A.4| ) implies that / = it; are the fixed points fot the 
first order flow of the quantity and as such they cannot be crossed in finite "time." 
The corresponding statement in the case of eq. ( 2.35| ) is prooven in Appendix (|C|). As an 
important consequence of this fact, the (true) energy of any solution interpolating between 
+v and —V is always positive, as we will see in the next subsection. 

Another useful form of £ is obtained by using the equations of motion ( |2.26| ) in ( 2.35| ) 
to put it in the form 



£ = ^^^^ mf - (v'u)?] 



ttR 



V{f) = 0. 



(2.42) 



2.2 The Bounce 

Next, we add a linear potential on the brane at y = ttR, of the form 

Vi = b{^-v), [b]=M^/'^. 
In this case, eq. ( |2.26| ) becomes: 

tan(7ri?a,) d^f - b = A / [f-v^] . 



(2.43) 



(2.44) 



The addition of ( 2.43| ) breaks the $ — > — $ symmetry and lifts the degeneracy between 
the two vacua, making one of them metastable. We will follow Coleman O and estimate 
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the decay rate of the metastable vacuum. We look for a "bounce" solution ^^{y,tE,x), 
i.e. a solution of the Euclidean field equations that interpolates between the true vacuum 



at small p = y t|, + and the false vacuum at large p. The bounce must have finite 
action relative to the false vacuum. Then the tunneling amplitude is given by: 

r = exp{-5s[$^] + 5s[$^]}. (2.45) 

Here S'^f'^] is the euclidean version of the action ( |2.1D evaluated on the field configuration 
<I>^ (the bounce) and (the false vacuum). 

We look for a bounce solution with 0(4)-symmetry, i.e. depending only on y and on the 
Euclidean radial coordinate p. The euclidean field equation in these coordinates reads: 

dl<^ + dl<i> + -dp<^ = (2.46) 

dy^\y^^ = m^^ -v^), (2.47) 
|^_^ = h. (2.48) 

and we look for a solution that approaches the true vacuum at /9 ~ and the false 
vacuum at p ~ oo. Following Coleman, we consider the symmetry breaking term as 
a perturbation: we approximate both the true and the false vacuum to be the same as 
the unperturbed ones (<I>(x,y) = ±t;), and moreover we set 6 = when solving the field 
equation. As a further approximation, we assume we are in the "thin wall" limit, in which 
we can neglect the last term in eq. ( p. 46 ). This is justified when the transition between 
the true and false vacuum takes place in a in a small region of width I}, around a radius 
po ^ h- Under these assumptions, the problem reduces to the one of the previous section, 
i.e. finding a domain wall solution centered around po- 



^^(P,y) 



-V < p < po, 

^kink{p - Po, y) P~Po (2-49) 
V P > Po, 



Requiring that the bounce has minimal action provides a variational problem for the 
parameter po. The bounce action is: 



^ = £ dyj dpp^\ [{d.^^f + {dy^^^f] - J dpp' (Fo(^>^(p,0)) - yi(^.^(p,7ri?))) . 

(2.50) 

Integrating by parts the bulk piece and using the field equations the above expression 
reduces to boundary terms. Approximating the solution as in ( 2.4S| ), we obtain: 



^ ^ -2bv^ + p^S^aii ■ (2.51) 



.3 



Here, S^aii is the energy stored in the wall. In our approximation can be thought as 
concentrated in a small region around po, and can be approximated by the total energy of 
the kink, i.e. the solution of eq. (2.26) : 

^ dx (V{^) - -CD— J ^^ = J^dx-{v'- fLki^)) . (2.52) 
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In the 4D regime, in which eqs. ( |2.29| ) and ( 2.30 ) hold, then after an integration by parts 



we obtain: 

f +00 



dx {dxfkinkf = -TTRfJ-v'^ = - VttRXv^ (2.53) 
Using this result in eq. ( 2.51| ) and minimizing the action with respect to po we find: 



p„ ^ ^ ^ (2-54) 

bv 



and the thin wall approximation holds if 

K m = — (2.55) 

which is the same condition [P finds in the purely 4D case, and that we would have obtained 
had we started from the 4D effective action in Section 2.1.1. 

From the previous discussion, it is clear that the thin- wall approximation gets worse 
and worse as we move away from the 4D regime, i.e. as RXv^ becomes large. In fact, as 
discussed earlier, from eqs. (^!3l|) and (|]3|) it follow that for RXv"^ >> 1 the width of the 



wall becomes much larger than the size of the region where the field profile has its largest 
variation (i.e. close to the / = 0), therefore the wall energy density gets spread over a 
larger and larger region. 

3. The warped case 

We are now going to repeat the steps in the previous section in a slice of AdS^ bounded 
by two branes at y = 0, ttR. We parametrize the metric as: 

ds^ = dy^ + e-'^^yr]^^dx^dx'' . (3.1) 

We will always assume a large warping, exp[A;7ri?] » 1. 
The field equation is 

- Akdy^ + e^'^f = . (3.2) 

with the same boundary conditions as before. The general solution for a mass eigenstate, 
□4^" = m?^, has the form: 



= e''y B2 [-e'y) (3.3) 

where B^{x) = a„iJu{x)+bmNiy[x) is an appropriate combination of Bessel functions, whose 
coefficients are to be determined along with the mass eigenvalues m, from the boundary 
conditions (^.3|-2.4). The correct linear combination is : 



and the resulting equation for the KK masses is: 

Bi{m/k) 2 



B2{m/k) 
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Typically the lowest KK mass is of order: 

mkk ~ A;e-^"^. (3.6) 

If we are in the 4D regime, when the lowest mode has mass mo << nikk-, expanding the 
Bessels in (|3^) we find: 



so the 4D regime demands that 



Xv^ <C yfce-^'^"^ . (3.8) 

In this 4d regime, the wave function of the zero mode is basically flat, whereas the 
^ — $ symmetry breaking term is redshifted by the warp factor b 6exp(— 4/c7ri?). 
The 4d Coleman expression for the bounce action is therefore valid and produces a huge 
enhancement of the lifetime of the false vacuum compared to the unwarped case. This is 
one of the main advantages in constructing metastable vacua in warped spaces. 

However, it would be very interesting to also understand the opposite regime, namely 
when /Uq = Af^ is much larger than the KK scale. Notice that in this regime the mass 
eigenstates are approximately given by the solutions of B2{m/k) = 0, since in this case the 



l.h.s. of eq. (3^) is large. 

Let us now look for a kink-like solution, depending on the coordinates x and y. From 
the flat case, we learned how to read-off an effective one-dimensional equation for the field 
at y = from the spectral equation. Repeating the argument that leads to ( |2.26| ), starting 
from the general bulk solution with correct boundary condition at y = ttR 



<^ix,y) = e^^y j dp ap 



rTt:^. — I ^(-' 0) ' (3-9) 

' Ox 

k 

k 



we obtain: 



'^2(f)-^^fc^iV2(^ 

Bi{m/k) 



B2{m/k) 



, d.f = Xf[f-v^] , (3.10) 
where /(x) = $(x,0). 

In the 4D limit we can take the first term in the expansion of the l.h.s. of eq. ( 3.1C| ), 
and we find again a second order equation, of the form: 

2kX 

{i-~7- 

which is the usual kink equation. 

Let us estimate the width of the kink in the opposite regime, Xv^ ^ ^tk/^- Using the 
same argument as in the previous section, and writing f{x) = ±v + lye^^'/'™ we find that 
Iw obeys: 

, B,{l/{kL)) _ , 



d'J= n _:l.nJ [f--'] (3-11) 
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which comparing with eq. (^]^) means that the maximal width is equal to the inverse mass 
of the lowest KK mode, i.e. of order ( |3.6| ). 

Let us now add the linear term ( |2.43| ) on the IR brane, as we did in the flat case. Making 
the same approximations as in Section 2.2 (treat Vi as a perturbation, and use the thin-wall 
approximation), we arrive at the following bounce action: 



[""dye-"'' f dpp'l [e''\d,<^y + {dy^^ 



^ = Jo '''' 

dpp' (yo(^^(p,0)) - e-4^-^6($^(p,7ri?) - v)) 

^ -2bv^e-'"'^'' + plS^aii, (3.13) 

where again we have assumed that ^^{p,y) = ^'^{p,y) = —v for p < pq, ^^{p,y) = 
^^{p,y) = +v for p > pq, and ^^{p,y) = ^kink for p « tq. Notice the appearance of the 
warp- factor in the first term of eq. ( |3.13| ) . S^aii is the same as in eq. ( |2.52| ), and it is 
localized on the brane at y = (we are neglecting the subleading contribution from the IR 
brane to the wall energy). 
Minimizing eq. ( 3.13| ) with respect to po we find: 



Swall AknR fo t a\ 

~ (3.14) 

and if Swaii is not too small this leads to an exponentially large radius of the vacuum 
bubble, and hence an exponentially small decay rate. 

We can give a crude estimate of Swaii as follows. Assume that fkink{x) can be approxi- 
mated piece- wise as: 

fkinki^) ~ ' ~lw < X — < 1^ (3.15) 

V X > Iw 

v 

Then evaluating the l.h.s. of ( |2.52| ) with this approximation we find 

S^all ~ 1^\V^ ^ ^^-^^^ 

In practice this may be an overestimate, since the linear regime assumed in ( p.l5| ) may not 
be valid for the whole width of the wall. But we can say that S^jaii is larger than just the 
contribution from the linear region: 

S-wall > Smin ~ hin^V^ (3-17) 



where lun is the region around the origin where the linear approximation ( |3.15 ) is justified 



It seems reasonable to believe that this region is independent of R for large enough R. In 
the flat case this region is of the order p^ = Xv^ for large pQ. We can repeat the same 
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analysis of Section (|2.1.2| ) to estimate the slope of the solution near x = 0. Assuming a 
behavior of the type f{x) ~ sin(x//o) we get the following equation: 



Mo (3.18) 



Now, let us assume that Mq ^ 1, so we can expand the Bessel functions evaluated in (klo)'^ 
(but not the ones evaluated in {klo)~^e'^^^ ). The quantity Ji{ie^'^^ /klQ)/Ni{ie^'^^ /Mq) 
is never small, since Ji{ix) has no zeros outside the origin. Using this fact, and expanding 
the Bessel functions of argument (A:Zo)~^ only, we obtain to lowest order: 



1 



2Xv^. (3.19) 



This result was obtained under the assumption Mq ^ 1, therefore it is valid if Af^ <^ k. 
This assumption is needed anyway, since we are in curved space, with a curvature scale of 
order k, and we are neglecting the backreaction of the scalar field <I> on the background, as 
well as the contribution of V{^) to the brane stress tensor. 
Using lo from (3.19) as an estimate of the width of the linear region, we get a a more 



conservative lower bound on S^aii from ( 3.17| ): 



S^^aii > loXv^ = y (3.20) 



up to 0(1) coefficients. If this is the case the size of the bounce is: 



Po > y exp[4fc7rfl] (3.21) 
and the thin-wall approximation holds if 

Po/lw : exp[3A;7ri?] > 1 , (3.22) 



which is easily satisfied due to the warp-factor. The decay rate is also exponentially sup- 
pressed: plugging ( 3.21 ) into (|3.13| ) we obtain: 



A2^;9 

Sb- -^exp[12kTTR] (3.23) 

which gives a huge lifetime r = e'^* even for moderate warping. This estimate shows that, 
in the warped case, we don't need to restrict to the 4-dimensional regime^ in order to have 
a small vacuum decay rate. 



®As this regime demands that Xv^ <^ k exp[—2knR], this would impose a very strong constraint on the 
model parameters. See however the next section for a different model. 
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There is an important omission in our previous discussion, the possible gravitational 
effects on the creation of the bubble. Indeed, Coleman and de Luccia showed |||] in the 4d 
context that gravitational effects are negligible only in the case 

^ « 1 , (3.24) 

where Aq is the radius such that the bubble radius equals the Schwarzchild radius. For 
the 4d version of our model it equals Aq = {IGGNbv/S)"^^"^ , where Gn = 1/Mp is the 
4d Newton constant. In our case and when the 4d approximation is valid, the two length 
scales scale with the warp factor as 

po exp[4A;7rii] po , Aq ^ exp[2A;7rii] Aq . (3.25) 

Then, neglecting factors of order one, gravity effects on the creation of the bubble are 
negligible when 

1 IXv^ 

exp[2kTTR] — \ — < 1 . (3.26) 
Alp \ 

If the 4d limit ( |3.8D is satisfied but ( |3.26| ) is violated, as shown in Q, there are two different 
cases. In the first, the metastable vacuum has positive energy whereas the true vacuum 
where we live has zero energy. Then gravity effects increase substantially the probability 
of tunneling. In the second case, the metastable vacuum has zero energy and tunnels 
into a negative energy stable vacuum. In this case, gravity effects increase the lifetime 
of the metastable vacuum. In the limit where /?o > 2Ao the bubble cannot form and 
the metastable vacuum becomes completely stable. This becomes therefore one important 
outcome of having a warped extra dimension, in the case where ( |3.25| ) is violated. 

4. Supersymmetric extension : the AdS-ISS model 

Recently, there was a renewed interest in metastable vacua from the point of view of 
supersymmetry breaking with further applications to gauge mediation models and 
moduli stabilization [13|. The proposal in P] used the electro- magnetic Seiberg duality to 



argue for the existence of metastable vacua in the supersymmetric QCD with a number of 
flavors A^c + 1 < < 3A^c/2. In the IR free magnetic description and before adding the 
effects of the (magnetic) gauge group, the model is described by the O'Raifeartaigh-type 
model 

W = hq^q - hjl'^Tr^ , (4.1) 

where q^ (q^) are the magnetic quarks (antiquarks), are the mesons, a,b = 1 ■ ■ ■ N are 
color indices and i,j = 1 ■ ■ ■ Nf are flavor ones. Supersymmetry is broken by the "rank 
condition", in the sense that the supersymmetry condition 

F$ = hqq - hjl^lNj , (4.2) 

where I^f is the A'^^ x Nf identity matrix, cannot be satisfied, since qq is a matrix of rank 
at most equal to A^ < Nf. One of the important requirements for the metastable vacuum 
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to be long-lived in the ISS model is e = jl/ A.m << 1, where is the Landau pole of the 
magnetic theory. From a string theory viewpoint one natural realization of the ISS 
model, in its magnetic description, is in terms of D3/D7 brane configurations, with the ISS 
gauge group realized on the D3 branes, with (anti) quarks coming from the D3-D7 sector 
and the magnetic mesons being the positions of a stack of D7 branes. 

The purpose of this section is to analyze in a field-theoretical example the effect that the 
warping of the internal space, generated by the branes, could have on the model. We model 
this effect by considering a five-dimensional supersymmetric model in a slice of AdS^ Q 
with the metric 

(4.3) 



dsl = e~^''^y^T]^udx^dx'' + dy^ 



with ISS gauge fields and the quarks, antiquarks confined to the UV boundary y = 
and the mesons promoted to a hypermultiplet (<I>i, ^2) propagating into the 5d bulk, with 
Z2 parities (-I-, — ). The mesons-quark coupling is localized on the UV brane, whereas we 
choose to put the linear term in the {Z2 even) mesons <I>i in the superpotential on the IR 
brane^. As we will show below, due to the exponential warp factor, the mass parameter 
fl will be redshifted such that the lifetime of the metastable vacuum becomes arbitrarily 
large. In a manifest 4d supersymmetric language |1C, 14], the Lagrangean describing the 
system is 

S = j d^xdy I j d^9 e-2*^2/($t$^ + ^t^^) + J d^Oe-^'^y {<^2dy<^i + h.c) 



+ 



j d^e {qU + qU) + j d^9 {hq^iq + Wnp{^i)+h.c) 
j d^e e^3fc^^( + h.c) 6iy - ttR) 



6{y) 



(4.4) 



where Wnp is the non-perturbative mesonic superpotential arising in the field direction 
where the mesons <l>i get vev's, give masses to the quarks (antiquarks) and generate the IR 
dynamics restoring supersymmetry. The (metastable) supersymmetry breaking becomes 
now a non-local effect and arises do to the impossibility, in the absence of Wnp, to solve 
the supersymmetric condition: 

e-^i^yF^, = -dy (e-^''y<^2) + (hqq + d^.Wnp) 6iy) - e'^'^-^ hfi^lNj5{y - ttR) . (4.5) 



In order to cancel the last term in eq. (|4.5| ), the Z2-odd mesons $2 acquires a non-trivial 
profile: 

^2 = e3'=(^--^)(/i/iV2)/iV,6(y). (4.6) 



If Wnp = supersymmetry is broken: with (|4.6| ), and using dye{y) = 2[5{y) — 5{y — ttR)], 



eq. ([4. 5]) becomes: 



e-^'^yp^, =5{y) 



hqq - e-^^^^ hfl^lNf 



(4.7) 



geometrical construction in a string context, similar in spirit, was proposed in 
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which cannot vanish due to the rank condition. Notice that the parameter which controls 
supersymmetry breaking is not /i, but rather 

^J'eff 6 A* ) Since QQ = e J-n ■ v^-°) 

The presence of Wnp restores supersymmetry by producing sources which do add up 
to zero. From the point of view of the bulk fields, in the metastable vacuum $i gets a 
boundary mass term 

^2 = /j2^^tg + q\q) = 2h^N e-^^'^^/i^ , (4.9) 

whereas in the supersymmetric vacuum it gets also localized nonperturbative interactions. 
The formally divergent terms 5(0) in ( |4.4| ) do not appear in physical quantities, as shown 
in various similar situations |jll|. 

Notice the close analogy of this model with the toy model analyzed in section 3: the 
symmetry breaking parameter is redshifted by a power of the scale factor. However in 
this model the validity of the 4D limit is automatic, and does nor require an additional 
fine tuning: the existence of a light mode for requires fi^ << A;exp[— 2/cr] , and from 
eq. ( |4.9D we see that this does not impose any strong constraint on h and /i, provided 
the warp factor is large. Therefore, since the 4D limit analysis holds, the smallness of 
the symmetry-breaking parameter due to the redshift leads immediately to an exponential 
enhancement of the lifetime of the metastable vacuum. There is one critical point to check: 
this conclusion is valid if the wave function of the lightest mode of $1 does not grow too 



fast in the IR and destroys the redshift of the mass term fl, transparent in (4^). In the 
limit where the 4d effective theory is valid, i.e. the lightest mode is much lighter than the 
KK masses m << ke~^'^^, its corresponding wave function reads approximatively 



2 

2kv 



+ d2 



2 

2ky 



4A;2 



(4.10) 



Boundary conditions determine then the mass spectrum to be given by the equation 

-2m^e-^'-^ = (^^ - ^) (4^ - ^e^^-^) . (4.11) 
Due to the validity of the 4D limit, we get the 4D result (see section 3) m? ~ 2A;/Xq and 



a corresponding wavefunction (4.10) which is constant in y to the leading order. In this 
case, the redshift of the mass parameter — > fj?e~^^'^^ is effective and produces a huge 
enhancement of the lifetime of the false vacuum. Notice that with respect to a 5d flat 
metric, the light mode is actually localized on the UV boundary. Since the KK modes and 
the linear term are localized on the IR boundary, this explains the enhancement of the 
lifetime of the metastable vacuum. 

Another interesting case, with the same matter content, is when the whole superpotential 
is localized on the UV boundary. In this case, there is no redshift of the mass parameter fi 
and generically no light mode. One way to obtain a light mode even in the case ke~'^^'^^ << 
« /c is to add a bulk mass for the hypermultiplet, which is tuned appropriately against 
the boundary mass. This is a tuning in a non-supersymmetric setup, but the tuning is 
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actually required and protected versus radiative corrections by supersymmetry |14, 16]. In 
this case the bulk mass mf, and the boundary masses fiQ, fi-j^ for the scalar component of 
$1, in the false vacuum, are given by 

'3 



2 



= - - c) , (4.12) 

where a = |c + 1/2|. Since we want to preserve in the first approximation the AdS^ 
geometry, we are interested in small backreaction of the scalar field and therefore small 
bulk mass q ~ 2. There is one interesting example of this type, with c = —5/2 and 
therefore zero bulk mass for <I>i, with non- vanishing brane localized masses. In this case 
we find a light scalar mode localized on the IR brane, with wave-function and mass given 
by 

$f^(x,y) ~ e-3'="^e^'=^(/.(x) , 



Qkh'^ {q'^q + q^) e-^^^^ . (4.13) 



The term exp(— S/cTri?), important in what follows, comes from normalization of the 4d 
kinetic term of the light mode (^{x). The four dimensional Lagrangean in this case is very 
close to the 4d ISS Lagrangean. Auxiliary fields are 

e-^'^yp^, = -dy{e-^^y^2) + {hqq - hfx^ + d^,Wnp)6{y) , 

Fg = e-^^^^(t)q , Fg = e-^'^^^qcf) . (4.14) 



Therefore, due to the wave- function in ( [4.13 ), the meson-quark coupling gets changed and 
become 

^-6knR^2^^^2 (|g|2 + |g~|2) _ (415) 

In the ISS vacuum, the quark vev's are as in 4d 

<i = f = { 1 , (4.16) 




Then ( [4.15 ) reproduces the light meson mode (4.13). In the SUSY vacuum in which 



mesons get vev's, quark masses are also redshifted by the same factor rrig = rv? 



exp(— 6A:7ri?)/i^|(/)p. Therefore the distance in field space between the ISS and the SUSY 
vacuum is greatly enhanced Ac/) = exp[3A;7ri?A^/(A'^y — N')\^^iss^ whereas the barrier re- 
mains unchanged T^eafc = -/V/^^/i^- Therefore the bounce action S\) in the triangular ap- 
proximation 5*5 ~ (A(/))^/VpeaA: fl^ and the lifetime of the false vacuum are accordingly 
increased 

12fc7rfliV 

5b ^ e ^/-^ 5b . (4.17) 
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However, as discussed in the previous section, a more detailed analysis of gravitational 
effects is needed in order to check if they are negligible. Again, if the metastable vacuum 
has zero energy whereas the stable vacuum has negative one, one expects the lifetime to 
be increased and eventually the false vacuum to become completely stable ||2[. 

Notice that for values /i^/i^ ~ k and by defining the mass scale on the IR brane with 
a dynamical scale A = /c exp(— /cvri?), we can rewrite qualitatively ( [4.13 ) in the suggestive 
way 

(4.18) 



where Mp is the 4d Planck mass. It is interesting to notice the analogy between ( 4.1^ ) 
and the scale of supersymmetry breaking in the observable sector in M = 1 supergravity 
with a gaugino condensation (AA) = A^ in a hidden sector, coupled gravitationally with 
the observable one. 

The models presented here can be interpreted from a holographic point of view. The 
metastable susy breaking can be understood in a purely four-dimensional way as arising 
from the infrared dynamics of a strongly coupled CFT sector, dual to the bulk geometry 
and the bulk fields <I>i,2- This CFT acts as a hidden sector, coupled to the quarks living 
on the UV brane. In the first example presented in this section (zero bulk mass), the light 
mode mediating the vacuum decay is localized on the UV brane, and from the point of 
view of the 4D theory it is a fundamental degree of freedom. The redshift of the mass 
parameter fl could be interpreted as the holographic version of the retrofitting discussed in 
|17, 19]. In the second example, in which the bulk field profile is given by eq. ( [4.13| ), the 
light mode is peaked on the IR brane and couples only gravitationally to the UV brane. 
In both cases, in the holographic 4D theory description the symmetry breaking occurs as 
an infrared effect, generating a hierarchy of scales like in eqs. (4^) and ( 4.18| ) 



In other types of models, in the nontrivial limit in which the boundary masses are large 
and the KK modes are expected to play a role in the bounce, there is no light mode anymore 
in the spectrum and the methods of Section 3 are needed in order to estimate the lifetime 
of the false vacuum. Finally, we would like to point out that there is nothing peculiar about 
the ISS model from the point of view of a phenomenological construction in a 5d warped 
space. Traditional O'Rafeartaigh models can be similarly discussed, with corresponding 
mass parameters and consequently scale of supersymmetry breaking redshifted to very 
small values. Since our main motivation was to understand the properties of the classical 
kink and bounce solutions, we refrain ourselves to discuss further here these applications. 
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A. Standard Kink solution 

Here, we remind the reader of the standard domain wall, or kink, solution. Consider a 
scalar field with quartic potential, 

Vm = ^{<S>'-v')' . (A.l) 



The one-dimensional field equation. 



q,"(x) = ^ = A$($2 _ v'^) , (A.2) 



with boundary conditions <!>(— oo) = —v, $(+oo) = v is solved by: 



<^kinkix) = vtanhfix, ^ = V ~- (^-^^ 
This is also a solution of the first order equation: 

S = ^ V{^) = . (A.4) 

This can be read as the conservation of energy equation of a point particle moving in the 
potential —V with vanishing total "energy" E. 
The total energy of the kink is 



E = 2 r dxV{^) = = —V{0)- . (A.5) 

J-oo 3 3 II 

B. Conserved energy 

In this appendix we derive the conserved "energy," eq. ( 2.35| ). We start from equation 
( 2.26 ) which we write in the form 

oo 

Y^and'^f-V'{f) = 0, (B.l) 

n=l 

where the a„ are defined by tan(7ri?x) = an^;^"""^. Next we multiply (B.l) by df and 
use the following identity 

2n-l 



dfd'^'^f = \dY^ (_i)P+iaP/52"-P/ (B.2) 

to get 



2 

p=i 



^ oo 2n— 1 

E «« E - V{f)] = 0. (B.3) 



n=l p=l 
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We deduce the conserved quantity 

oo 2n— 1 

2 



^ oo 2n— 1 



(B.4) 



n=l p=l 

Symbolically the sum '^'^^^(^—l^P+^d^fd'^"'~^f can be written as 



9l + ^2 

where di = dx^- The first term in ( p.4D can thus be put in the form 
d d 

-(tan(7ri?5i) + iSin.{-KRd2))f{xi)f{x2)\xi=x 



2=X 



di + d2 
Now we use 

tan(7ri?(9i) + tan(7ri?92) = [1 — tan(7ri?9i) tan(7ri?92)] tan(7ri?(9i + 52)) 
and {di + f {xi) f {x2)\x^=x2=x = <9"/^, which gives 

tan(7ri?(ai + 82)) f {xi) f {x2)\xi=x2=x = t&n{7rRd)f. 
Collecting all the terms we get the final expression 

tan{-KRd) 



ttRS = - 
2 



d 



[idff-idtani7rRd)ff]-Vif). 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B. 



C. The 5D Kink near the extrema of the scalar potential 

Here we analyze the behavior of the kink in flat space, close to the extrema of the potential, 
/ = 0, zizv. In particular we show that a solution with zero "energy" £, i.e. satisfying eq. 
(|2.42 ), cannot cross from a region where |/| < v to another one where |/| > v. 



We have already shown in Section 2.1.2 that when / ~ f the solution to eq. ( p. 261) is 
exponential, 

f ^v±rieM±^/L], ^ tan(^] = 2Xv'' , (C.l) 



where is a constant. One can check that the above ansatz satisfies the condition £ = 
to lowest order in r]: inserting (|C.lJ ) in ( ^.35 ) and keeping terms quadratic in rj we obtain 
(for any choice of signs in ( |C.lD ): 



1 tan(7ri?(9) 
2 

/ 1 tan(27ri?// 



t an (vri?//^ ) -r— e 



If, 



±2x/L 

ttR 



2 2TTR/k 



{l - tan^ {it R/l^)) 



MP 
ttR 



2+2x/l^ 



>r] e 



— tan — - 2Xv^ )■ — 
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where in the last hne we used the identity: 



tan 2z 



2tanz 
1 — tan^ z 



(C.2) 



Each solution of the type ( C.l ) approaches ibw as |x| — > oo. One can ask whether it is 
possible for a solution to approach (and cross) |/| = u at a finite value x = xq. A priori, 
one can take : 

f v + r] sinh[(2; - xo)/lw] x ^ xq , (C.3) 

as a solution to the linearized kink equation ( ^.261) with the desired property to cross f = v 
at X = xq. However, let us compute the conserved energy for ( p.3|) : 



£ 



1 tan(7ri?a) 

2 irRd 



— cosh[(x - xo)/k 



(tan(7ri?9)9sinh[(x — xo)/lw]Y 



Xv rj sinh [(x — xo)//«,]. 



(C.4) 



We will use the following formal identity: for any differential operator 0{d) constructed 
with a function 0{k) which has an expansion containing only even powers of k, (such as 
the two operators appearing in the above expression) we have: 

0{d) sinh kx = 0{k) sinh kx, 0{d) cosh kx = 0{k) cosh kx . (C-5) 

Using this fact, and some manipulation of the hyperbolic functions, we arrive at : 



£ 



1 

4E 



2 (t^R 

1 + tan^ 



^ 2ttr] 



(C.6) 



therefore a solution that crosses / = if cannot have a zero value of £. 

On the contrary, a zero energy solution can cross / = at some finite value of x. Close 
to / = the solution of the linearized equation has now the form (see eq. ( |2.32| )): 



f{x) ~ r/sin[(x — xq)/Io], tanh[7ri?//o] = Af Iq, (C.7) 
Inserting this in eq. ( |B.9D and performing the same steps that led to eq. ( |C.6| ) we obtain: 



1 



4/2 



£ = i—r, [1- tanh^ — 



ttR 



^ 27rR} 



1] 



AttR' 



(C.8) 



The last term is ??- independent, and comes from the non-zero value of V{f) at / = 0. For 
an appropriate choice of r/, we can make £ vanish: 



Xv^ 



1 



1 — tanh^ 



ttR 



Xv^ \ 

^ 2TTRj 



£ = 



(C.9) 



Notice that this argument does not require r] to be small: the validity of the linearized 
approximation made in eq. ( p.Tj) holds for arbitrary r], as long as x is close enough to xq. 
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